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We have measured the size effect in nonlocal Gilbert relaxation rate in FM(t fm) / 
Cu (5nm) [/ Pt (2nm)] / Al(2nm) heterostructures, FM = { Ni 8 iFei 9 , Co 60 Fe 2 oB 2 o, 
pure Co}. Common behavior is observed for three FM layers, where the additional 
relaxation obeys both a strict inverse power law dependence AG = Kt n , n = — 1.04± 
0.06 and a similar magnitude K = 224 ± 40 Mhz • nm. As the tested FM layers span 
an order of magnitude in spin diffusion length Xsdl, the results are in support of 
spin diffusion, rather than nonlocal resistivity, as the origin of the effect. 
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The primary materials parameter which describes the temporal response of magnetization 
M to applied fields H is the Gilbert damping parameter a, or relaxation rate G = \j\M s a. 
Understanding of the Gilbert relaxation, particularly in structures of reduced dimension, is 
an essential question for optimizing the high speed / Ghz response of nanoscale magnetic 
devices. 

Experiments over the last decade have established that the Gilbert relaxation of ferro- 
magnetic ultrathin films exhibits a size effect, some component of which is nonlocal. Both 
ot(tFNi) = ao + a (tpM) and G{tpM) = Gq + G {Ifm) increase severalfold with decreasing FM 
film thickness tpM, from near-bulk values ao,Go for tpM ~ 20 nm. Moreover, the damp- 
ing size effect can have a nonlocal contribution responsive to layers or scattering centers 
removed, through a nonmagnetic (NM) layer, from the precessing FM. Contributed Gilbert 
relaxation has been seen from other FM layers^ as well as from heavy-element scattering 
layers such as Pt.- 

The nonlocal damping size effect is strongly reminiscent of the electrical resistivity in 
ferromagnetic ultrathin films. Electrical resistivity p is size-dependent by a similar factor 
over a similar range oitpM] the resistivity p(£fm) is similarly nonlocal, dependent upon layers 
not in direct contact It is prima facie plausible that the nonlocal damping and nonlocal 
electrical resistivity share a common origin in momentum scattering (with relaxation time 
t m) by overlayers. If the nonlocal damping arises from nonlocal scattering r^ 1 , however, 
there should be a marked dependence upon FM layer type. Damping in materials with 
short spin diffusion length Xsdl is thought to be proportional to (ref.-); the claim for 
" resistivity- like" damping has been made explicitly for Ni 81 Feig by Ingvarsson^ et al. For FM 
with a long Xsdl, on the other hand, relaxation G is either nearly constant with temperature 
or "conductivity- like," scaling as tm- 

Interpretation of the nonlocal damping size effect has centered instead on a spin current 
model^ advanced by Tserkovnyak et air-. An explicit prediction of this model is that the 
magnitude of the nonlocal Gilbert relaxation rate AG is only weakly dependent upon the 
FM layer type. The effect has been calculated^ as 



, where the effective spin mixing conductance gljj/S is given in units of channels per area. 
Ab-initio calculations predict a very weak materials dependence for the interfacial parameters 
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g^/S, with ±10% difference in systems as different as Fe/Au and Co/Cu, and negligible 
dependence on interfacial mixing.— 

Individual measurements exist of the spin mixing conductance, through the damping, 
in FM systems NigiFeig^, CcM, and CoFeB^. However, these experiments do not share 
a common methodology, which makes a numerical comparison of the results problematic, 
especially given that Gilbert damping estimates are to some extent model-dependent^. In 
our experiments, we have taken care to isolate the nonlocal damping contribution due to Pt 
overlayers only, controlling for growth effects, interfacial intermixing, and inhomogeneous 
losses. The only variable in our comparison of nonlocal damping AG(tpM), to the extent 
possible, has been the identity of the FM layer. 

Gilbert damping a has been measured through ferromagnetic resonance (FMR) from 
a;/27r = 2-24 Ghz using a broadband coplanar waveguide (CPW) with broad center conduc- 
tor width w = 400/xm, using field modulation and lock-in detection of the transmitted signal 
to enhance sensitivity. The Gilbert damping has been separated from inhomogeneous broad- 
ening in the films measured using the well-known relation AH pp (uj) = AH + (2/ y/3j au/\^\. 
We have fit spectra to Lorenzian derivatives with Dysonian components at each frequency, 
for each film, to extract the linewidth AH PP and resonance field H res ; a has been extracted 
using linear fits to AH(lo). 

For the films, six series of heterostructures were deposited of the form Si/ Si02/ 
X/ FM(W)/ Cu(3nm)[/Pt(3nm)]/ Al(3nm), FM = { Ni 81 Fe 19 ("Py"), Co 60 Fe 20 B 2 o 
("CoFeB"), pure Co}, and t FM = 2.5, 3.5, 6.0, 10.0, 17.5, 30.0 nm, for 36 heterostruc- 
tures included in the study. For each ferromagnetic layer type FM, one thickness series tpM 
was deposited with the Pt overlayer and one thickness series tpM was deposited without the 
Pt overlayer. This makes it possible to record the additional damping Aa(tpM) introduced 
by the Pt overlayer alone, independent of size effects present in the FM/Cu layers deposited 
below. In the case of pure Co, a X=Ta(5nm)/Cu(5nm) underlayer was necessary to sta- 
bilize low-linewidth films, otherwise, depositions were carried out directly upon the in-situ 
ion-cleaned substrate. 

Field-for-resonance data are presented in Figured) The main panel shows u(H B ) data 
for NigiFeig(tFM)- Note that there is a size effect in oo(H^): the thinner films have a 
substantially lower resonance frequency. For tpM = 2.5 nm, the resonance frequency is 
depressed by ~5 Ghz from ~20 Ghz resonance Hg — 4 kOe. The behavior is fitted through 
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the Kittel relation (lines) u(H%) = \j\J[H 1 ^ + H K j [AttM!^ + H% + H K j, and the inset 
shows a summary of extracted AttMI^^fm) data for the three different FM layers. Samples 
with (open symbols) and without (closed symbols) Pt overlayers show negligible differences. 
Linear fits according to AitM^f (tpM) — AirM s — (2K S /M S ) t FM allow the extraction of bulk 
magnetization AnM s and surface anisotropy K s ] we find AnM Pv = 10.7 kG, AnM^ oFeB = 
11.8 kG, 47rMf° = 18.3 kG, and Kfv = 0.69 erg/cm 2 , K? oFeB = 0.69 erg/cm 2 , = 
1.04 erg/cm 2 . The value of g L /2 = \^\/(e/mc), \y\ = 2tt ■ (2.799 Mhz/Oe) • (g L /2) is found 
from the Kittel fits subject to this choice, yielding g Py = 2.09, gL° FeB — 2.07, g£° = 2.15. 
The AirM s and values, taken to be size- independent, are in good agreement with bulk 
values. 

FMR linewidth as a function of frequency AH pp (uj) is plotted in Figure El The data 
for Py show a near-proportionality, with negligble inhomogeneous component AHq < 4 Oe 
even for the the thinnest layers, facilitating the extraction of intrinsic damping parameter 
a. The size effect in in aitphi) accounts for an increase by a factor of ~ 3, from a Py = 
0.0067 (G Py = 105 Mhz) for the thickest films (t FM = 30.0 nm) to a = 0.021 for the 
thinnest films (tFM = 2.5 nm). The inset shows the line shapes for films with and without 
Pt, illustrating the broadening without significant frequency shift or significant change in 
peak asymmetry. 

A similar analysis has been carried through for CoFeB and Co (not pictured). Larger 
inhomogeneous linewidths are observed for pure Co, but homogeneous linewidth still ex- 
ceeds inhomogeneous linewidth by a factor of three over the frequency range studied, and 
inhomogeneous linewidths agree within experimental error for the thinnest films with and 
without Pt overlayers. We extract for these films a$ oFeB = 0.0065 (Gfi oFeB = 111 Mhz) 
and olq° = 0.0085 (Gq° = 234 Mhz). The latter value is in very good agreement with the 
average of easy- and hard-axis values for epitaxial FCC Co films measured up to 90 Ghz, 
G%° = 225 Mhz.i£ 

We isolate the effect of Pt overlayers on the damping size effect in Figure [3] Values 
of a have been fitted for each deposited heterostructure: each FM type, at each £_fm, 
for films with and without Pt overlayers. We take the difference Aa(tpM) for identical 
FM(t^jv/)/Cu(5nm)/Al(3nm) depositions with and without the insertion of Pt(3nm) after 
the Cu deposition. Data, as shown on the logarithmic plot in the main panel, are found 
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to obey a power law Aa(tFM) = Kt n , with n = -1.04 ± 0.06. This is excellent agreement 
with an inverse thickness dependence Aa(t F M) = KpM/tFM, where the prefactor clearly 
depends on the FM layer, highest for Py and lowest for Co. Note that efforts to extract 
Aa(tFAi) = Kt n without the FM(tpA/)/Cu baselines would meet with significant errors; 
numerical fits to a(t FM ) = KtpM n for the FM(t^ A ./)/Cu/Pt structures yield exponents 
n ~ 1.4. 

Expressing now the additional Gilbert relaxation as AG(t F M) = \~f\M s Aa(t F M) = 
\l FM \Mf M K FM /t FM , we plot AG • t FM in Figure H We find AG ■ t Py = 192 ± 40 Mhz, 
AG • tcoFeB = 265 ± 40 Mhz, and AG • t Co = 216 ± 40 Mhz. The similarity of values for 
AG ■ tpM is in good agreement with predictions of the spin pumping model in Equation [H 
given that interfacial spin mixing parameters are nearly equal in different systems. 

The similarity of the AG • t F M values for the different FM layers is, however, at odds 
with expectations from the "resistivity-like" mechanism. In Figure HI inset, we show the 
dependence of AG -tpu upon the tabulated Xsdl of these layers from Ref^. It can be seen 
that Xsdl * s roughly an order of magnitude longer than it is for the other two FM layers, 
Py and CoFeB, but the contribution of Pt overlayers to damping is very close to their 
average. Since under the resistivity mechanism, only Py and CoFeB should be susceptible 
to a resistivity contribution in Aa(t F M), the results imply that the contribution of Pt to 
the nonlocal damping size effect has a separate origin. 

Finally, we compare the magnitude of the nonlocal damping size effect with that pre- 
dicted by the spin pumping model in Ref.— . According to AG • t F M = {'y^h/An = 
25.69 Mhz • nm 3 ((7i/2) 2 (gljf/S^j, our experimental AG • t F M and gi data yield effective 
spin mixing conductances gtjf/S [Py/Cu/Pt] = 6.8 nm~ 2 , gljf/S [Co/Cu/Pt] = 7.3 nm" 2 , 
and gtjf/S [CoFeB/Cu/Pt] = 9.6 nm . The Sharvin-corrected form, in the realistic limit 
of \% DL » ijyli is (gtjf/S)-' = (gf /N /S)-' - foftg/S)- 1 + 2e 2 ^V t N + (g% l/N JS)-\ 
Using conductances 14.1nm~ 2 (Co/Cu), 15.0nm~ 2 (Cu), 211nm -2 (bulk pc u , tN = 3nm), 35 
nm" 2 (Cu/Pt) would predict a theoretical gtjfth /S [Co/Cu/Pt] = 14.1 nm -2 . Reconciling 
theory and experiment would require an order of magnitude larger pc u — 20yuf2 • cm, likely 
not physical. 

To summarize, a common methodology, controlling for damping size effects and intermix- 
ing in single films, has allowed us to compare the nonlocal damping size effect in different 
FM layers. We observe, for Cu/Pt overlayers, the same power law in thickness £-l-04±0.06^ 
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the same materials independence, but roughly half the magnitude that predicted by the spin 
pumping theory of Tserkovnyak—. The rough independence on FM spin diffusion length, 
shown here for the first time, argues against a resistivity-based interpretation for the effect. 
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FIG. 1. Fields for resonance u(H B ) for in-plane FMR, FM=Ni 8 iFei 9 , 2.5 nm < t FM < 30.0 nm; 
solid lines are Kittel fits. Inset: 4-7T for all three FM/Cu, with and without Pt overlayers. 




FIG. 2. Frequency-dependent peak-to-peak FMR linewidth AH pp (lo) for FM=NigiFei9, tFM 
as noted, films with Pt overlayers. Inset: lineshapes and fits for films with and without Pt, 
FM=Ni 8 iFei 9 , CoFeB. 
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FIG. 3. Inset: a no p^(^fm) and ctp^ for Py, after linear fits to data in Figure [2j Main panel: 
Aa(tpM) = apt(^FAf) — « no pt(^FM) for Py, CoFeB, and Co. The slopes express the power law 
exponent n = -1.04 ± 0.06. 
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FIG. 4. The additional nonlocal relaxation due to Pt overlayers, expressed as a Gilbert relaxation 
rate - thickness product AG • tFM for Py, CoFeB, and Co. Inset: dependence of AG • tpM ° n sphi 
diffusion length Xsdl as tabulated mil. 
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